Rayleigh-Benard Convection with a Radial Ramp in Plate Separation 



M. R. PauH and M. C. Cross 
Department of Physics, California Institute of Technology 114-36, Pasadena, California 91125 

P. F. Fischer 

Mathematics and Computer Science Division, Argonne National Laboratory, Argonne, Illinois 60439 

(Dated: February 8, 2008) 

Pattern formation in Rayleigh-Benard convection in a large-aspect-ratio cylinder with a radial 
ramp in the plate separation is studied analytically and numerically by performing numerical sim- 
ulations of the Boussinesq equations. A horizontal mean flow and a vertical large scale counterflow 
are quantified and used to understand the pattern wavenumber. Our results suggest that the mean 
flow, generated by amplitude gradients, plays an important role in the roll compression observed 
as the control parameter is increased. Near threshold the mean flow has a quadrupole dependence 
with a single vortex in each quadrant while away from threshold the mean flow exhibits an octupole 
dependence with a counter-rotating pair of vortices in each quadrant. This is confirmed analytically 
using the amplitude equation and Cross-Newell mean fiow equation. By performing numerical ex- 
periments the large scale counterfiow is also found to aid in the roll compression away from threshold 
but to a much lesser degree. Our results yield an understanding of the pattern wavenumbers ob- 
served in experiment away from threshold and suggest that near threshold the mean flow and large 
scale counterfiow are not responsible for the observed shift to smaller than critical wavenumbers. 



I. INTRODUCTION 

Rayleigh-Benard convection in a thin horizontal fluid 
layer heated from below is a canonical example of pat- 
tern formation in a continuous dissipative system far from 
equilibrium §. Under various conditions, wavenumbcr 
selection mechanisms have been identified ||, ^, |4| that 
reduce the band of stable wavenumbers, sometimes to a 
single value. One such selection mechanism occurs when 
there is a one-dimensional spatial variation or ramping 
of the control parameter e; where e = {R — Rc)/Rc, 
R is the Rayleigh number, and Rc its critical value 
^ H, 0, H, ^, 0. This can be accomplished by varying 
the plate separation d such that e goes from e = eo > 
in the bulk of the layer (i.e., the unramped region) to 
e < as a lateral boundary is approached. It is expected 
that in the idealized case of an infinitely gradual one- 
dimensional ramp the wavenumber will equal kc at the 
position where the layer depth yields Rc- It has been 
shown under very general conditions, that as long as the 
layer becomes critical somewhere along the ramp this is 
sufficient to fix the wavenumber fc^ in the bulk and over 
the rest of the ramp Q . For slightly supercritical condi- 
tions it is expected that the selected wavenumber in the 
bulk can be expressed as 
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(1) 



where kc = kc = 3.117 |9j and a depends on the Prandtl 
number, a, and the specifics of the ramp. 

Recent Rayleigh-Benard convection experiments |l^ 
in a cylindrical cell with a two-dimensional radial ramp 
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in plate separation have generated intriguing results. 
Specifically, the plate separation as a function of radius 
used in experiment is, using the layer depth d to nondi- 
mensionalize. 
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where Sr is a constant and the radius values tq and ri 
are the locations where the ramp begins and ends, re- 
spectively. The ramp always extends to the sidewalk A 
schematic of the cosine ramp is shown in Fig. |^; note 
that To and ri are geometric constants but rc, the loca- 
tion where the plate separation yields Rc, is a function 
of the ramp shape and eo such that, for a ramp given by 
Eq. d), drc/deo > 0. 
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FIG. 1: A vertical cross-section of a cylindrical convec- 
tion layer with a radial ramp in plate separation, ro defines 
where the ramp begins, ri defines where the ramp ends and 
Tc is where the plate separation corresponds to the critical 
Rayleigh number. The ramp shown is a cosine ramp given 
by Eq. (j^. For presentation purposes we show a steep ramp 
with Sr ~ 0.25 (this domain is not used in the simulations). 

Experiments using the cosine ramp defined by Eq. iMj 
have yielded unexpected results for the wavenumber |||, 
pd| . Mean pattern wavenumber measurements (using the 
Fourier methods discussed in |l^) yielded kc = 2.97 < 
kc- Additionally, measurements of the local wavenumber 
defined at each position in space (method discussed in 
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|l3[) displayed interesting variation as eo is increased. 
For the time independent patterns the bulk region, r < 
To, contains approximately straight parallel rolls. Near 
threshold, eo < 0.048, a centered egg-shaped domain of 
convection rolls with small wavenumber, k < kc, extends 
through the convection cell with the long-axis parallel to 
the roll axes. A dramatic roll expansion from fc « 3.6 
at the edge where the ramp begins to fc w 2.6 in the 
center of the domain is observed. As eo increases, the 
wavenumber field evolves into a domain characterized by 
large wavenumbers extending through the layer with the 
long- axis perpendicular to the roll axes (see Fig. 3 of [^). 
Similar experiments without a ramp do not exhibit this 
wavenumber behavior |l^, |l6|. 

In the experiments with a radial ramp in plate sepa- 
ration, for eo ^ 0.03, time dependent states were found 
through the repeated formation of defects via an Eckhaus 
mechanism consistent with the local wavenumbers cross- 
ing the Eckhaus stability boundary for an ideal infinite 
layer of two-dimensional rolls. Also in the ramped exper- 
iments, for eo ^ 0.18, defects were formed via a skewed 
varicose mechanism consistent with the local wavenum- 
ber exceeding the skewed varicose stability boundary for 
an ideal infinite layer of two-dimensional rolls similar 
to what has been observed in experiments in unramped 
cylindrical domains with rigid sidewalls [ p^ . 

It has been suggested that these features may be the 
result of the interaction of the convective roll pattern 
and weak large scale flows The visualization and 
quantification of these large scale flows is not possible 
in the current generation of experiments. However, we 
are able to make these measurements by performing full 
numerical simulations with a new spectral element code 
(discussed further in Section III). We utilize the com- 
plete knowledge of the flow fleld together with analytical 
results valid near threshold to explore this further. 



II. LARGE SCALE FLOWS 

In this work the terminology large scale flows is used 
to describe flows that extend over distances larger than 
that of the convection roll scale. We would like to distin- 
guish between two different large scale flows: large scale 
counterflow and mean flow. 



A. Large Scale Counterflow 

In the presence of a spatial ramp in plate separation 
a large scale counterflow is present for all values of the 
bulk control parameter eo, including eo < 0. Warm fluid 
ascends the ramp eventually reaching the sidewall and 
is forced to flow back toward the center of the domain 
over the cold top wall causing it to descend resulting 
in a large zone of circulation in the vertical plane over 
the ramp. The magnitude of the large scale counterflow 
depends upon the specifics of the ramp and is roughly 



independent of cq and a [jisj . Figure || illustrates this 
with a vertical slice from a three-dimensional numerical 
simulation where the entire convection layer is subcriti- 
cal. As shown for this subcritical case the fluid motion of 
the large scale counterflow generates axisymmetric con- 
vection near the base of the ramp which extends a couple 
of roll widths toward the center of the domain. For the 
gradual ramp used in experiment the large scale counter- 
flow is small in magnitude and has not been measured. 



FIG. 2: Velocity vectors illustrating the large scale counter- 
flow for a vertical slice of a cylindrical convection layer for 
subcritical conditions, eo — —0.063. The ramp parameters 
are ro — 7.66, ri = 10, and Sr — 0.15. A steep ramp is shown 
to clearly illustrate the flow. Solid vertical lines indicate the 
boundaries of the spectral element grid used in the simulation, 
only a portion of the layer is shown emphasizing the ramped 
region. 



B. Mean Flow 

In discussing the mean flow it will be convenient to first 
present the governing heat and fluid equations. The ve- 
locity u, temperature T, and pressure p, evolve according 
to the Boussinesq equations, 

a^^ (dt + u-V^u = -\'p + RTz + S/^u, (3) 
(dt + U'V'jT = V^T, (4) 
V.M = 0, (5) 



where dt indicates time differentiation, and z is a unit 
vector in the vertical direction opposite to gravity. The 
equations are nondimensionalized in the standard man- 
ner using the layer depth d, the vertical diffusion time for 
heat Ty = (P' I K where k is the thermal diffusivity, and AT 
the temperature difference between the top and bottom 
surfaces, as the length, time, and temperature scales, re- 
spectively. The lower and upper surfaces are no-slip and 
are held at constant temperature. The sidewalls are also 
no-slip and, unless otherwise noted, are insulating. 

The mean flow fleld, U{x,y), is the horizontal veloc- 
ity integrated over the depth and originates from the 
Reynolds stress induced by pattern distortions. As il- 
lustrated by the fluid equations, Eqs. (||) and (||), it is 
evident that the pressure is not an independent dynamic 
variable. The pressure is determined implicitly to enforce 
incompressibility. 



V^p=-a-^W' 



Rd.T. 



(6) 
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Focussing on the nonlinear Reynolds stress term and 
rewriting the pressure as p = Po{x, y) + p{x, y, z) yields, 



Po{x,y)^a-^ / rfxWln(l/|r-/|)(V'. 



(7) 

where {■) ^ represents an average in the z direction. In 
Eq. (0) the ln(l/|r - r'\) is not exact, in order to be 
more precise the finite system Green's function would 
be required; however the long range behavior persists. 
This gives a contribution to the pressure that depends on 
distant parts of the convection pattern. The Poiseuille- 
like flow driven by this pressure field subtracts from the 
Reynolds stress induced flow leading to a divergence free 
horizontal flow that can be described in terms of a verti- 
cal vorticity. 

Near threshold an explicit expression for the mean flow 

is ra 



U{x, y) = --fkV^. (^^2) - V^poix, y) 



(8) 



where 7 is a coupling constant given by 7 = 0.42cr~^((T -I- 
0.34)(cr -t- 0.51)~^, is the convection amplitude nor- 
malized so that the convective heat flow per unit area 
relative to the conducted heat flow at Rc is \A\'^R/Rc, Po 



and Vj^ is the 
for the complete 



is a slowly varying pressure, see Ec 
horizontal gradient operator [see 
analysis and more details]. The mean flow is important 
not because of its strength; under most conditions the 
magnitude of the mean flow is substantially smaller than 
the magnitude of the roll flow making it extremely diffi- 
cult to quantify experimentally. The mean flow is impor- 
tant because it is a nonlocal effect acting over large dis- 
tances (many roll widths) and changes important general 
predictions of the phase equation |]l9j . The mean flow is 
driven by roll curvature, roll compression and gradients 
in the convection amplitude. The resulting mean flow 
advects the pattern giving an additional slow time de- 
pendence. It is important to note, that unlike the long 
range counterflow, the magnitude of the mean flow van- 
ishes when the convection layer becomes critical, \U\ ~ eo 
for eo < 1. 



III. NUMERICAL SIMULATION 

We have performed full numerical simulations of the 
governing fluid and heat equations, Eqs. (|1)-(|1), in a 
cylindrical geometry with a radial ramp in plate sepa- 
ration using a parallel spectral element algorithm (de- 
scribed in detail elsewhere see Q for related 
applications). 

For discussion purposes it will be convenient to define 
Cartesian {x, y) and polar (r, 6) coordinates centered on 
a mid-depth horizontal slice of a cylindrical convection 
layer containing a field of straight parallel x-rolls with 
wavevector k = koX as shown in Fig. H. The x-axis is 



perpendicular to the roll axes, the y-axis is parallel to 
the roll axes and measures the angle from the positive 
X-axis. 




FIG. 3: Cartesian {x,y) and polar (r,9) coordinates defined 
on a mid-depth horizontal cross-section of a cylindrical con- 
vection layer containing a field of x-rolls described by fe = koX. 

We have investigated the results found in experiment 
by performing simulations on a ramped cylindrical con- 
vection layer for a variety of scenarios and initial con- 
ditions. Simulations were performed over the range 
eo ^ 0.2 and for simulation times oi tf > r/j where 



Th 



— ^2 



is the time required for heat to diffuse horizon- 
tally across the bulk region of the layer which has been 
suggested as the earliest time scale for the flow field to 
reach equilibrium . 

The mean flow present in the simulation flow fields, 
Us{x,y), is investigated by calculating the depth aver- 
aged horizontal velocity. 



Us{x,y) = / u±dz 



(9) 



where u± is the horizontal velocity field. Furthermore it 
will be convenient to work with the vorticity potential, 
C, defined as 



Vic 



z. V I X U. 



(10) 



where uJz is the vertical vorticity and is the horizontal 
Laplacian. 



IV. ANALYTICAL DEVELOPMENT 

Near threshold, assuming straight parallel rolls, it is 
possible to approximately determine coz analytically. It 
will be convenient to start from uj^ given by the vertical 
component of the curl of Eq. (||) , 



(fc|AP 



(11) 
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Consider a cylindrical convection layer with a radial ramp 
in plate separation containing a field of x-roUs given by 
k = koX. The amplitude can be represented for large eo, 
using an adiabatic approximation, as = e{r)/go for 
e > and \A\'^ = for e(r) < 0, making the amplitude a 
function of radius only \A\'^ = /('')• This approximation 
is good except for the kink at Tc where e = 0. Inserting 
|Ap — f{r) into Eq. ( [lT| ) yields, after some manipulation, 
the following expression for the vertical vorticity. 



£\A\^ _ l d\A['' 



r dr 



sin 20. 



(12) 



To correct for nonadiabaticity and to smooth |A(r)|^ 
near Tc, the one-dimensional time independent amplitude 
equation is solved, 

O = e(r)A + eoW0^-.go|A|2A, (13) 
where angular derivatives have been assumed small, 



= 0.148, go = 0.6995 
e(r) is determined by 



0.0047cr- 



0.0083(7"^ and 
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eo, 



eo{h^-Kl)/{l~hl), 



r < ro 
r > ro 



(14) 



where he = hire)- In Eq. ( [Til ) is dominated by radial 
derivatives so Eq. dl2) is still a good approximation for 
the case = f{r^, d/dr d/dr. Due to the angular 
dependence of the nonadiabaticity \A\'^ is now 9 depen- 
dent which will induce higher angular harmonics in Uz- 
We will neglect these higher harmonics, assume a sin 20 
dependence and approximately evaluate the magnitude 
using Eqs. (|l|) and (|l|) at 6* = 7r/4. Equation (^ can 
be rewritten in a more convenient form as, 



d'^A 

Qj:2 



A^ 



(15) 



where 



— -1/2 

— (eo ' 



Cocos0)f and A = {cQ/gof^A |1 



(since the amplitude goes to zero it can be shown that 
A is now a real quantity (2^). Equation ( |l5| ) is solved 
numerically using the boundary conditions dfA = at 
f = 0, and A = at f = fi . 

These analytical results are now used to investigate the 
vertical vorticity generation in a large radially ramped 
cylindrical convection layer. We start by looking at the 
configuration used in the recent experiments with input 
parameters: tq = 42.29, ri ~ 101.33, 5r = 0.036 and 
cr = 0.87. 

When eo^/^Co~^(f'c — "Tq) < I the amplitude A^{r) is 
unable to adiabatically follow the ramp, this nonadia- 
baticity results in a considerable deviation from e(r)/go 
as shown in Fig. 0. However, when eo^/^^o~^(?'c~?'o) ^ 1 
the amplitude ^ (r) follows e{r)/go adiabatically almost 
over the entire ramp except for the small kink at Vc as 
shown in Fig. ^. The structure of ujz depends upon this 
adiabaticity and is shown for various values of eo in Fig. ^ 



where the sin 29 dependence has been removed by choos- 
ing 9 = 7r/4. 

If the 9 dependence is included it is evident from Fig. ^ 
that the vertical vorticity has a quadrupole angular struc- 
ture for small eo, i.e. four lobes of alternating positive 
and negative vorticity with one lobe per quadrant, and 
makes a transition to an octupole angular dependence for 
larger eo, octupole in the sense of an inner, r < Tc, and 
outer, r > Tc, quadrupole. In addition, since drc/deo > 
there is a radial shift of the vorticity curves as eo is in- 
creased. 

In all cases the amplitude A'^{r) decreases monoton- 
ically with r and as a result — r^^ d\A\'^/dr > thus 
generating only positive vorticity. However, the term 
d'^\A\'^ /dr'^ can be of either sign and is responsible for 
the quadrupole and octupole angular structure in ujz- 
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FIG. 4: The solution of Eq. ^ plotted as A^(r) for ro = 
42.29, rc = 47.56, n = 101.33, 5r = 0.036, a = 0.87 and 
eo = 4.20 X 10~^. Also shown for comparison is e{r)/ go- 
As eo approaches zero the nonadiabaticity of A^(r) in- 
creases until A^ (r) exhibits a quadratic fall-off with r for 
r ^ ro resulting in Wz(r) > for all r. 

The mean flow generated by these vorticity distribu- 
tions is determined by solving Eq. ( |ll| ) with the bound- 
ary condition C('''i) ~ 0. The vorticity potential is re- 
lated to the mean flow in polar coordinates by {Ur, Ug) = 
{r~^dgC, —drC)- The vorticity potential is expanded ra- 
dially in second order Bessel functions while maintain- 
ing the sin 20 angular dependence. Of particular interest 
is the mean flow perpendicular to the convection rolls, 
Ur{9 = 0) or equivalently Ux{y — 0), which is shown in 
Fig. 0. 

As expected, regions of negative and positive vorticity 
yield corresponding negative and positive values of the 
mean flow. As eo vanishes U > for all r providing a 
mechanism for roll expansion in the bulk. For larger eo 
the mean flow becomes larger in magnitude and increas- 
ingly negative for r < ro providing a mechanism for roll 
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FIG. 7: Variation of the mean flow, Ux{y = 0) with eo as 
FIG. 5: The solution of Eq. (|l|) plotted as A'^{r) for ro = determined from Eq. @. Shown exphcitly are curves for 
42.29, rc = 80.05, n = 101.33, 5r = 0.036, a = 0.87 and = 7.2 x 10"^, 4.2 x 10-^ and 1.3 x 10"^ labelled (a)-(c), 

eo = 0.17. Also shown for comparison is e{r)/go. respectively with the parameters ro = 42.29, n = 101.33, 

Sr = 0.036 and a = 0.87. 
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5,0, = eoT-\ and r-i - 19.65a(a + 0.5117)-i [|. 
Assuming that for eg ^ 1 the wavenumber is approxi- 
mately k Ki kc everywhere and that the rolls are exposed 
to a constant mean flow the wavenumber change over the 
bulk can be expressed as 



Afc = k{ro) - k{r = 0) = UkcDu^Q. 



(17) 



For example, for curve (c) in Fig. M the maximum value 
of the mean flow is U — 4.92 x 10 which yields a small 
roll expansion of Ak = 0.0035. If the mean flow were 
solely responsible for the dramatic roll expansion seen in 
experiment of Afc « 1.0 (eo = 0.012 see Fig. 3(a) of §) a 
mean flow of C/ w 0.014 would be required, which is not 
found in the analytic results. 



FIG. 6: Dependence of the radial variation of cuz = fr l 
with eo as determined analytically from Eqs. ( [l^ ) and ( [l^ ) 
illustrating the evolution from a mean flow with a quadrupole 
dependence to an octupole dependence as eo increases. Shown 
explicitly are luz curves for eo = 7.2 x 10"^ 4.2 x 10'^ and 
1.3 X 10"^ with the parameters ro = 42.29, n = 101.33, 
Sr ~ 0.036 and a = 0.87 labelled (a)-(c), respectively. 

compression. 

To make the connection between mean flow and 
wavenumber quantitative it is noted that the wavenum- 
ber variation resulting from a mean flow across a field 
on x-roUs can be determined from the one-dimensional 
phase equation, 

Ud^cj) ^ D^\d^,(b (16) 
where the wavenumber is the gradient of the phase, k = 



V. DISCUSSION 

The large scale flows discussed in Section || cannot be 
measured in current experiments placing us in a unique 
position to use the complete flow field information from 
our full numerical simulations of Eqs . (|^ )-(|^) together 
with the analytical results of Section p^V| to investigate 
how the mean fiow and the large scale counterfiow in- 
duce wavenumber distortions and the variation of this 
distortion with Rayleigh number. 

It is computationally expensive to perform full three- 
dimensional numerical simulations for the very large sys- 
tem used in experiment. We have, however, performed 
a variety of simulations for radially ramped cylindrical 
convection layers. The full three-dimensional simulations 
are of smaller spatial extent with the precise ramp defined 
by Eq. (^ and the specific input parameters: rg = 11.31, 
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ri = 20, 5r = 0.036, and a = 0.87. Two-dimensional sim- 
ulations of a vertical slice of a three-dimensional domain 
(see Fig. |l|) were also performed for both the large ex- 
perimental configuration and the smaller computational 
domain just described. Three-dimensional simulations 
were also conducted without a large scale counterflow by 
a specific choice of ramp parameters that will be dis- 
cussed below. 

Initially we consider prescribed x-roll initial conditions 
given by fc = k^x. Other initial conditions such as ran- 
dom thermal perturbations or initial x-roUs of varying 
wavenumbers were also investigated and found not to af- 
fect the final pattern wavenumber or any of the conclu- 
sions drawn. Simulations were performed for eg = 0.025, 
0.054, 0.113, and 0.171. Figure ^ compares the wavenum- 
bers found in these simulations with recent experiments 
and will be discussed in detail below. 
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subcritical region as the supercriticality of the bulk in- 
creases. All simulations settled to a time independent 
state. 
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FIG. 8: Comparison of the mean wavenumber variation as a 
function of the bulk Rayleigh number (i.e. for r < ro), -Ro, 
between simulation (solid lines with symbols) and experiment 
(dashed-dotted line) |^]. Unless otherwise noted ro = 11.31, 
ri = 20, 6r = 0.036, and a — 0.87. The symbols represent; 
(o) three-dimensional simulations, (□) two-dimensional simu- 
lations, (>) three-dimensional simulations for a specific ramp 
construction without a large scale counterflow, and (o) two- 
dimensional simulations with ro = 42.29, and ri — 101.33. 
Dark solid lines denote the approximate location of the neu- 
tral (N) and Eckhaus (E) stability boundaries for an ideal 
infinite layer of parallel rolls. 

The final patterns in the simulations maintain the 
x-roll configuration imposed by the initial conditions. 
Figure ^ displays the final pattern observed for three- 
dimensional simulations with ep — 0.025 in panel (a) and 
eo = 0.171 in panel (b). Figure ^a illustrates that near 
threshold the convection rolls exhibit very little curva- 
ture indicating that the assumption of straight parallel 
x-rolls in Section [V is valid. There is more roll curva- 
ture apparent in Fig. ^ as would be expected for larger 
eo- Figure 1^ also illustrates the decreasing size of the 




FIG. 9: Final convection patterns for eo = 0.025 and 0.17 are 
shown in panels (a) and (b) , respectively. Shaded contours of 
the thermal perturbation are shown with dark regions repre- 
senting cool descending fluid and light regions warm ascend- 
ing fluid. The inner dotted circle indicates where the ramp 
begins, ro, and the outer dotted circle indicates where the 
convection layer becomes critical, r^. Simulation parameters, 
ro = 11.31, ri = 20, 5r = 0.036, a = 0.87. 

It is illustrative to compare the analytical results of 
Section fv| with the results of simulation. Figure [lo|a 



displays A^(r) for the case eo ~ 0.025, as determined by 
Eq. (y_5|). A significant nonadiabaticity is present for this 
case as shown by the deviation of A'^{r) from e{r)/go- 
For the ramped domain used in simulation the distance 
^c(eo) ~ fo is smaller than in the larger domain with a 
more shallow ramp used in experiment. This results in 
the presence of more nonadiabaticity in the simulations 
when compared to experimental results at the same con- 
trol parameter. This is beneficial because this allows 
the exploration of highly nonadiabatic situations without 
having to perform the task of simulating near the convec- 
tive threshold, which becomes computationally difficult 
because of the diverging time scales. 

A comparison between theory and simulation of the 
vertical vorticity and the resulting mean fiow is shown 
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in Fig. liol b and c. The theoretical predictions are based 
on the amphtude variation caused when straight parallel 
convection rolls encounter a radial ramp in plate sepa- 
ration as discussed in Section [V. For both the vertical 
vorticity and the mean flow the comparison is made in 
the absence of any adjustable parameters. For the verti- 
cal vorticity calculated in simulation an angular average, 
weighted by sin 20, is used for the comparison. The agree- 
ment between theory and simulation is quite good. This 
illustrates quantitatively that the major source of verti- 
cal vorticity and mean flow is indeed the variation in the 
convective amplitude caused by the radial ramp in plate 
separation. Over the bulk of the domain the mean flow is 
negative and very small in magnitude with a maximum 
value of U^y = 0) = -7.0 x lO"* and by Eq. ^ the 
wavenumber variation would be extremely small in agree- 
ment with the near constant bulk wavenumbers found in 
simulation. 

A similar comparison between theory and simulation 
is made in Fig. ^ for eg = 0.171. As shown in Fig. |Tl|a 
A'^(r) is much more able to follow the ramp, e(r)/(/o, and 
exhibits very little nonadiabaticity except for the kink 
near Vc- This results in a much stronger negative vertical 
vorticity in the bulk which in turn yields a larger negative 
mean flow as shown if Fig. pT| b and c. The agreement be- 
tween theory and simulation for the vertical vorticity is 
still quite good. The discrepancy in the mean flow com- 
parison may be due to the fact that as eo increases other 
mean flow sources such as roll curvature, see Fig. 
become important. 

Figure |l^ illustrates the octupole structure in the vor- 
ticity potential in panel (a) and the roll compression 
occurring in the bulk by plotting contours of the local 
wavenumber in panel (b) for cq = 0.171. As illustrated in 
panel (a) the mean flow has signiflcant structure over the 
ramped region as well as extending into the subcritical 
region of the layer, r > Tc- It has also been suggested that 
the mean flow extends into a subcritical region in related 
experiments implementing "flnned" boundaries |2^ . 

The vorticity potential displays an octupole structure 
containing a pair of counter-rotating vortices in each 
quadrant. The inner quadrupole is localized around tq 
where gradients in the amplitude of convection occur as 
the ramp in plate separation begins. The direction of 
rotation of the inner quadrupole causes a focussing of 
the mean flow into the bulk region of the domain and 
is responsible for the larger wavenumbers found as eo is 
increased as shown by the (o) curve in Fig. |^. 

To make the connection between mean flow and 
wavenumber quantitative Eq. (|l^) is applied to the sim- 
ulation results in the form. 



(18) 



Figure |13| a illustrates the wavenumber variation, fc(r) for 
r < tq, found in simulation by simply measuring the dis- 
tance between roll boundaries and makes evident the roll 
compression, k{r = 0) > k{ro). Figure [l3| b compares the 
mean flow calculated from simulation with the predicted 
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FIG. 10: Panel (a) shows the solution of Eq. (nSf) plotted as 
yl^(r), shown for comparison is e(r)/go. Panel lb) compares 
the vertical vorticity found analytically from Eq. (y_5|) with an 
angular average, weighted by sin 26, of the vertical vorticity 
from simulation. Panel (c) compares the mean flow found 
analytically from Eq. (^|) with the mean flow from simulation 
flowing along the x-axis &t y = Q. Parameters are ro = 11.31, 
rc = 13.20, n = 20.0, 5r = 0.036, a = 0.87 and eo = 0.025. 



value of the mean flow required to produce the wavenum- 
ber variation shown in Fig. [Tsj a using Eq. (|l|). The agree- 
ment is good and the discrepancy near ro, which is con- 
tained within one roll wavelength from where the ramp 
begins, is expected because the influence of the ramp was 
not included in Eq. (|l6|). This illustrates quantitatively 
that the mean flow compresses the rolls in the bulk of 
the domain. 

As mentioned earlier, the mean flow vanishes as eg ap- 
proaches critical whereas the large scale counterflow is 
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FIG. 11: Panel (a) shows the solution ol Eq. ( p^ ) plotted 
as A^ir), shown for comparison is e{r)/go. Panel b compares 
the vertical vorticity found analytically from Eq. (^5|) with an 
angular average, weighted by sin 26, of the vertical vorticity 
from simulation. Panel (c) compares the mean flow found 
analytically from Eq. (^|) with the mean flow from simulation 
flowing along the x-axis at y = 0. Parameters are ro — 11.31, 
rc = 13.20, n = 20.0, Sr = 0.036, a = 0.87 and eo = 0.171. 



present for all eo and therefore could play a role near 
threshold in the determination of the final convection 
pattern. In order to gain further insight into this pos- 
sibility a radial ramp was constructed that did not drive 
a large scale counterflow. This was accomplished by set- 
ting the temperature of the ramped surface, T'b(r), to 
the value of the linear conduction profile at that height, 
Tbir) — h{r). This ramp, therefore, does not bend the 
isotherms which is the source of the large scale coun- 
terflow. The wavenumber variation for these Simula- 



I 



0.06 



0.00 



-0 .06 



3.3 



3.2 



3.1 




FIG. 12: Contours of the vorticity potential (, panel (a), 
(light indicates counterclockwise rotation drawn with solid 
contours and dark indicates clockwise rotation drawn with 
dashed contours) and the corresponding local wavenumber 
distribution, k, panel (b). The magnitude of the mean flow 
is approximately 2% of the magnitude of the velocity fleld, 
|[/ij|/|M| w 0.02. Local wavenumber distributions are shown 
only in the bulk, r < vq. The inner dotted circle indicates 
where the ramp begins, ro, and the outer dotted circle indi- 
cates where the convection layer becomes critical, Vc- Simu- 
lation parameters, ro — 11.31, ri = 20, 6r = 0.036, a = 0.87 
and eo = 0.171 {Rq = 2000). 



tions, see curve labelled with (>) in Fig. does not dif- 
fer strongly from the simulations with a ramp produc- 
ing large scale counterflow, see the curve labelled with 
(o). The similarity in wavenumber results is strongest 
for small eq suggesting that the large scale counterflow is 
not responsible for the shift of the critical wavenumber 
to smaller values as seen in experiment. 

To study the large scale counterflow further two- 
dimensional simulations were also performed, corre- 
sponding to a vertical slice of the domain considered 
thus far, and in addition to a more spatially extended 
domain as used in experiment. In two dimensions the 
mean flow is absent however the large scale counterflow 
persists. As shown by the (o) and (□) curves in Fig. ^ 
the wavenumbers measured in the two-dimensional sim- 
ulations are not compressed to the same extent as eo in- 
creases as in the three-dimensional simulations with both 
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FIG. 13: Panel (a), the variation in the local wavenumber 
along the positive x-axis, or equivalently k{r) at = 0. 
Panel (b), a comparison of the mean flow from simulation 
(solid line) with the predicted value calculated from Eq. ( p^ ) 
using the wavenumber variation from panel (a). Simulation 
parameters, ro = 11.31, ri = 20, 5r = 0.036, a = 0.87 and 
eo = 0.171 (Ro = 2000). 

mean flow and large scale counterflow present. As ex- 
pected, the wavenumbers found in the two-dimensional 
simulations are also independent of aspect ratio. Addi- 
tionally, for small eq the wavenumber found in simulation 
does not deviate markedly from its critical value suggest- 
ing that the large scale counterflow is not responsible 
for the wavenumber shift observed in experiment near 
threshold regardless of the spatial extent of the domain. 

We also investigated the possibility that the vertical 
large scale counterflow could bifurcate into a horizontal 
flow similar to the mean flow in the presence of a slight 
spatial asymmetry. This was accomplished by giving the 
ramped domain used in the full three-dimensional simu- 
lations an eccentricity of e « 0.8 for a variety of ramps 
0.036 < 6r < 0.25 and simulating over a range of subcrit- 
ical and supercritical conditions. For all of the scenarios 
tested the large scale counterflow remained vertical and 
did not undergo any significant changes. 

Lastly, the possibility of wavenumber pinning was 
studied by varying the aspect ratio in increments of less 
than half of a roll width for both the two- and three- 
dimensional domains. In all of the scenarios tested the 



final pattern wavenumbers were not appreciably affected 
by these small changes in aspect ratio. 



VI. CONCLUSION 

We have analytically and numerically investigated pat- 
tern formation in a cylindrical convection layer with a 
radial ramp in plate separation. In particular, we have 
studied quantitatively the effects of two large scale fiows; 
large scale counterflow and mean flow. These large scale 
flows are important theoretically yet are extremely diffi- 
cult to measure experimentally. 

Our results suggest that the mean flow plays an im- 
portant role in the observed pattern wavenumber and is 
generated in a novel way, by the spatial variation of \A\'^ 
driven by the variation of e rather than the more usual 
variations in roll curvature and wavenumber. The mean 
flow sources are quantifled analytically and agreement is 
found with numerical results. The geometric structure 
and magnitude of the mean flow is used to explain quan- 
titatively the wavenumber variation found in the simula- 
tions. 

The large scale counterflow is investigated numerically 
and our results indicate a small roll compression effect 
away from threshold. In particular, although the large 
scale counterflow is present at and near threshold it 
does not appear responsible for the dramatic wavenum- 
ber shift to values less than critical as seen in experiment. 

Although it is too expensive computationally to sim- 
ulate the very large systems used in the experiments 
we can use our quantitative understanding of the ramp- 
generated mean flow, validated by the simulations at 
smaller aspect ratio, to extrapolate our results to these 
larger systems. Furthermore our two-dimensional sim- 
ulations in sizes equal to the experimental ones allow 
us to estimate the effect of the large scale counterflow 
on the wavenumber distribution. Despite these exhaus- 
tive efforts, we are unable to reproduce the large shift 
to smaller wavenumbers observed near threshold in the 
experiments, and the physical origin of these results re- 
mains a mystery. 
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